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Abstract 

Bose statistics imply a substantial enhancement at small angles for light 

scattering off a cold, Bose gas. The enhancement increases dramatically at the 

Bose-Einstein temperature. This phenomenon could be utilized to eliminate 

almost entirely the heating of the gas by a weak probe light beam. 
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A salient feature of Bose statistics is that the transition probability to a particular final 
state is proportional to n/ + 1, where rif is the number of bosons hitherto in that final state. 
For photons, this is the origin of stimulated emission due to the presence of an external 
electromagnetic field. The analog for the motion of bosonic atoms is the enhancement of 
scattering into already occupied states. This is the essential ingredient in the dynamics of 
the formation of a Bose-Einstein condensate as the temperature drops below the critical 
value, T c . A boson transition probability is also linear in the number of identical particles in 
the initial state. The present note points out how these factors effect simple photon-atom 
scattering in the Born approximation (e.g. for weak light). At low temperatures, there is 
substantial enhancement for small photon scattering angles (corresponding to small momen- 
tum transfers) from stimulated scattering of atoms into already occupied low momentum 
states. Below T c , this enhancement dramatically shifts to even smaller angles due to the 
contribution of transitions from occupied atomic momentum p ^ states to the Bose con- 
densate p = state and vice versa. Aside from seeing Bose statistics in action, this may 
have a practical implication for the optical investigation of cold atoms. The incident light 
intensity — and consequent total heating of the gas — can be reduced by several orders of 
magnitude for a given, desired signal strength by putting the detector at very small angle. 
This may allow non-destructive measurements on super-cold gases. 

I treat here the simplest case of a uniform medium and study the photon scattering rate 
per unit time and volume into a given solid angle. A further simplifying approximation is 
that the atoms are treated as an ideal, i.e. non-interacting, Bose gas, so that the multiparticle 
states are described by occupancies of single particle states. Whether appropriate in detail to 
a particular atom trap depends on trap shape and size, atomic density, and beam wavelength 
and size. However, the same qualitative features will emerge [[![]. 

The underlying process is just Rayleigh scattering (or perhaps on resonance to enhance 
the rate). I focus on the angular dependence of the rate for fixed incident wavelength. 
With very sub-Doppler temperatures T (and with a very narrow light frequency band if 
using a resonance — narrow compared to the atomic line width) , there is no relevant energy 
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dependence in the scattering amplitude even though we sum over possible initial and final 
atomic kinetic energies. This is because the range of energies involved in the thermal sum 
over atoms is negligible compared to the energy scale of variation of the scattering amplitude. 
The only angular dependence in the underlying "unstimulated" rate comes from the photon 
polarizations: cos 2 6 for polarizations in the scattering plane, 1 for perpendicular, |(1 + 
cos 2 9) for unpolarized - - all of which I treat as 1 for small angles. (The polarization 
difference may have some potential interest because it vanishes as 6 — > and, hence, weights 
angles differently.) 

For a given, pure, many-particle quantum state with n« atoms in an initial one-particle 
state and rif atoms initially in the one-particle state into which one of the rii is scattered, 
the scattering rate is proportional to n^nf + 1). The thermally averaged rate is, therefore, 
proportional to (rij(n/ + 1)), where the brackets denote the thermal expectation. For i ^ f, 
which is guaranteed by having 0^0, this is equal to (rii)((nf) + 1) if we ignore interactions 
between the atoms themselves. Hence, the needed thermal information is simply the thermal 
average occupancies of the single particle states i and /. It is simplest to evaluate these 
in terms of a chemical potential, which is then adjusted to produce the desired mean total 
number of particles. Note, however, that for bosons at very low temperatures, if any state 
has a large occupancy, the expected fluctuations in the occupancy are comparably large. 
The actual, physical situation of interest will determine which description (fixed particle 
number or fixed chemical potential) is more convenient. 

It is useful to divide the scattering rate into three terms, (a) The total "unstimulated" 
rate, i.e., the 1 term in rif + 1, is normalized to N totat ■ 1, where N tota i is the total target 
atom density and where the Rayleigh rate is normalized to I. (b) The scattering from 
p ^ occupied states to p ^ occupied states is represented by a double integral over 
two Bose thermal occupation factors (-1). And (c) below T c , there is scattering from p ^ 
occupied states to the p = state, with occupation N (t), and vice versa; these two pieces 
give identical contributions. Because the integrals over p's converge as p — > 0, they do not 
account for a macroscopic occupancy of the p = state, which must be treated separately. 
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Of course, there is no N Q ■ N Q term at non-zero angle. Note that the "unstimulated" rate is 
linear in the density, while the stimulated terms are quadratic. 

Choose units such that H = c = 1. The convenient third choice of units is mc 2 /csT c = 1, 
where m is the atom mass. Let A be the dimensionless measure of the momentum transfer 
to the photon, i.e. expressed in the natural units of the problem. Then 

. ... sin fa/ n , . 

where is the photon scattering angle, k is the magnitude of its momentum in these units, 
and v is its frequency. For laser light scattering off alkali atoms condensing at temperatures 
of order 10~ 7 K, 9 is of order 0.1 to 0.01 times A. (If the conversion factor between and A 
were much smaller, e.g. with yet lower atom density and, hence, lower T c , the phenomena 
discussed below might be limited to inaccessibly small angles.) 

Define a scaled temperature r = T/T c . Then i?(A,r), the total scattering rate normal- 
ized to the Rayleigh rate off a density of N totah is given by 

R(A,t) = 1 (2a) 

+N -i f ^VW(p-p'-a) 

total J [ A -l( r ) eP a/2r - l][A- 1 (r)ef' 2 /2r _ 1] ^ °> 
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N tota J [\-\r)e^-l\ 



+2 



where 



N total = J d* P (e^ - I)- 1 = \{2,f'\(\) 
N (t) = (1 - r 3 / 2 )iV totai e(l - r) 



A(r) = e^ T),T 

= 1 " { lL [C( ^ )]2(r " 1)2 + ° i{T ' 1)3)}0(T " 1} 

and where /x(r) is the chemical potential and the Riemann C(3/2) ~ 2.61. The coefficient 
in A(r) of (r — 1) above r = 1 is determined by requiring N total (which in these units is 



the critical density) to be independent of r. The virtue of the dimensionless variables A, r, 
and N tota i is that a single function R(A,t) describes all physical situations, and interesting 
phenomena occur when the variables are roughly 0{\). The actual light frequency, atom 
density, and atom mass then determine the translation of (A, r) into (0,T). 

The integral in term (2b) cannot be evaluated in closed form. After using the 5-function 
and integrating over angles, one is faced with a one-dimensional integral, depending on 
A and r, that must be performed numerically. Before addressing the results, note that, 
despite the 6-functions in N a and A, both R(A,r) and dR/dr are continuous functions of 
r. (The latter is not immediately obvious but can be confirmed by an explicit evaluation 
of the contribution to dR/dr from the r-dependence of A(r) at r = 1; this can be done 
analytically. The discontinuity in this contribution to the r derivative is exactly cancelled 
by the r-dependence of term (2c). The r-dependence arising from the e p2 / r 's is clearly 
continuous.) 

The Bose final state enhancement of the total scattering, i.e., integrated over all angles, 
is a factor of N tota i, independent of r. What changes with r is the angular dependence of that 
enhancement. As r decreases, it is concentrated at smaller and smaller angles — because 
the thermal distribution is at lower and lower momenta. What happens as we go below 
r = 1 is that the finite fraction of the Bose enhancement involving the p = Bose-Einstein 
condensate has a different angular distribution, concentrated at somewhat smaller angles. It 
is the singular nature of the Bose distribution as p — > and r — > 1 + that keeps R and dR/dr 
formally continuous. Nevertheless, the numerical evaluation reveals a dramatic crossover for 
sufficiently small A. 

The most interesting ranges of parameters are displayed in figure (1). Above the critical 
temperature, A's of order 1 or below show a Bose enhancement (solid lines) of a factor of two 
to five. Below the critical temperature, for A's like 0.1 or below, the enhancement can jump 
to several hundred or thousand with only a ten to twenty percent drop in temperature. The 
dashed curves are the contribution to the normalized rate coming from scatterings involving 
the p = states, either as targets or final states. These scatterings account for a substantial 
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fraction but by no means all of the enhancement. At fixed A, as T — > 0, the enhancement 
must ultimately disappear. (It is then concentrated into yet smaller angles). 

The above effects are really just simple kinematics that follow from the Bose occupation 
numbers. However, they have a potential practical application in cold atom optics. By 
collecting scattered light at small angles, one could significantly reduce the incident intensity 
needed to produce a given signal strength. The thermal average of the enhanced small angle 
scatterings do not heat the gas at all. They collectively remove precisely as much energy 
as they add because they involve all possible transitions from the thermal distribution back 
into the thermal distribution. The isotropic component (2a) is always present, and it does 
always contribute to heating. However, if the rate at some milliradians is 500 times greater 
than at 90 degrees, then one could reduce the incident intensity and, hence, total heating 
rate by 500 by moving the detector to small angles and still have the same signal strength. 

The present discussion makes the simplifying assumption that the target gas is uniform 
over the beam size. One practical effect of finite sample size is that diffraction limits the 
smallest meaningful angles. For current experiments 0, the conversion of A to a laboratory 
scattering angle puts the interesting angles at the border of what might be observed. A 
small atom trap may also produce significant variation in the position space support of the 
wave functions. This can be used, of course, to advantage in identifying the ground state 
occupancy 0. The formulae presented here are more appropriate to a trap that is large 
compared to the beam size or more like an ideal box. Nevertheless, the behavior for small, 
harmonic traps should be similar Jl]. There are again three classes of states: unoccupied 
high energy states, a thermal occupation distribution of low energy excited states, and, 
finally, the ground state, with a significant occupancy. And the "unstimulated" rate is 
linear in the density while the stimulated rate is quadratic. One salient difference between 
trapped atoms and a uniform medium is that, with an external potential, momentum is not 
conserved In an infinite, uniform medium, the super-stimulated ground-state-to-ground- 
state transitions correspond to zero photon scattering angle. The corresponding transitions 
in a trap are analogs of the Mossbauer effect and are dominant for photon momentum 
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transfers up to about double the largest typical ground state atom momentum. It is only 
beyond the corresponding maximum scattering angle that trap ground-state-to-excited-state 
and excited-state-to-excited-state scattering are the largest effects. Note, however, that if 
the light scattering is performed after a significantly long time interval of free expansion, 
i.e. after suddenly turning off the trapping fields 0, the Bose-stimulated enhancement 
will be confined to increasinly smaller angles. This is because, under free propagation, 
different momentum components become spatially separated. So the overlap between the 
scattered atom, of a particular final momentum, and previously existing atoms with that 
same momentum decreases rapidly with time. 

Many recent theoretical discussions of optical effects in low temperature Bose gases have 
focused on far subtler issues, such as coherence in forward scattering at low |3j or high 
densities (relative to the optical wavelength) or with pulsed light 0, or the effects of finite 
trap size |lj or transition line shapes [0]. (For a comprehensive review, see ||.) Although the 
fact that stimulated scattering occurs at small angles has often been mentioned in passing, 
a quantitative estimate of the magnitudes, angles, and sharpness with temperature of the 
effect had not been given. As emphasized in this context by Javanainen |J , to manifest any 
dramatic consequence of Bose statistics, the situation must involve a substantial number of 
identical particles that can be in the same place. For the present discussion to be relevant, 
the figure of merit is the number of atoms in the volume defined by the inverse of the 
momentum transfer. 
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FIGURES 

FIG. 1. The scattered light, total enhancement factor, R(A,T/T C ), including Bose stimulated 
atomic motion (solid lines) versus temperature at different reduced angles A (see eqn. (1)). The 
dashed curves are the condensate contributions (i.e., from eq. (2c)). 
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